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In this paper, we study the piezoresistive response and its anisotropy for a bond percolation 
model of granular metals. Both effective medium results and numerical Monte Carlo calculations of 
finite simple cubic networks show that the piezoresistive anisotropy is a strongly dependent function 
of bond probability p and of bond conductance distribution width Ag. We find that piezoresistive 
anisotropy is strongly suppressed as p is reduced and/or Ag is enhanced and that it vanishes at the 
percolation thresold p = pc. We argue that a measurement of the piezoresistive anisotropy could be 
a sensitive tool to estimate critical metallic concentrations in real granular metals. 
PACS numbers: 72.20.Fr, 72.60.+g, 72.80.Ng 



I. INTRODUCTION 



Granular metals are composite materials where a conducting phase made of metallic particles or clusters with mean 
sizes ranging from about 10 A to 200 A or more is randomly dispersed into an insulating glassy material. Depending 
on the nature of the insulating and conducting phases, their relative concentrations and the typical size of the metallic 
granules, these naaterials show a quite rich phenomenology and peculiar properties such as metal-insulator transitions, 
giant Hall effect,EI and variable-range hopping type of transport at low temperatures for insulating samples.EI 

The main ingredients governing most of the granular metal properties are the concentration x of metal in the 
composite, the activation energy Ec and the tunneling parameter K = 2d/^, where d is the tunneling distance and 
^ is the localization length. By varing x, the compiosite undergoes a metal-insulator transition at a critical value Xc 
which can be interpeted as the percolation thresold.u For x > Xc, connected metallic grains form a macroscopic cluster 
which extends over the whole sample and current can flow easily from one end to another when a voltage difference is 
applied. In this regime, transport has a metallic character identified by a positive temperature coefficient of resistivity 
(TCR). When x < Xc, transport is dominated by tunneling and the resistivity increases as the temperature is lowered 
(so TCR is negative). In this regime, the activation energy Ec plays a fundamental role since it defines the scale of 
energy an electron should overcome to successfully tunnel from one grain to another. It is this parameter that can 
give rise, together with K , to variable-range hopping type of transport at sufficiently low temperatures.Ba 

The experimental determination of Xc is however not as clear cut as the above discussion seems to suggest. For 
granular metals with low values of Ec and x close to the percolative thresold, TCR can change sign as the temperature 
is varied in such a way that, TCR < for T < T* and TCR > for T > T*, where T* is some temperature 
characteristic of the materialB As pointed out in Ref. H, the nonmonotonic behavior of TCR makes the determination 
of Xc ambiguous and alternative methods should be developed. 

In this paper we study the piezoresistive response of granular metals and show how its anisotropy is strictly related 
to the vicinity to the percolation critical point. Piezoresistance, i.e., the variation of resistance R or conductance G 
upon an applied strain field e, is an effect common also in bulk metals. If we consider an isotropic conducting material 
having bulk resistivity p, length I, and cross-sectional area A, then the resistance R is given by: 

The fractional change of R, SR/R, due to a strain e = 51/1 is given by: 

SR _ dp SI SA 

If v is the Poisson's ratio of the conductor, then it is possible to show that SA/A — ~2ve so that: 

GF= ^ + (l + 2z.), (3) 

ep 

where we have defined the total piezoresistive gauge factor GF = (5i?/ei?.0i It is clear that the total piezoresistive 
response is made of an intrinsic, Spjep, and a geometric, 1 -I- 2v, component. For bulk metals, GF is about 2 ami, 
since v is usually between 0.2 and 0.4, the geometrical effect is in this case more important than the intrinsic one.Q 
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A different situation is encountered instead in the so-called thick-film resistors, a particular class of granular metals 
typically made of au02, Ir02, or Bi2Ru2 07 granules embedded in a glassy matrix, for which GF values as large as ~ 30 
has been reported,El so that the intrinsic piezoresistive effect is much larger than the geometrical one. Such high strain- 
sensitivity resistances ha"\ie|-been exploited to manufacture piezoresistive sensor devices successfully used for pressure 
and force measurements. Lro Due to their commercial applications, thick-film resistors are the granular materials 
for which the piezoresistive efjEftet has been studied best and it is thought that intcr-grain tunneling processes are 
responsible for their GF values .lH13 To illustrate this point, consider for example the conductance g of a simple tunneling 
process between two metallic grains: g ~ exp(— X). Under an applied strain, the tunneling distance d is changed to 
d + 5d = d{l -I- e) where e = 5d/dJif the strain. This modifies the tunneling parameter, K K ^ 6K = K{1 -f e) and 
consequently the conductance g.Ej The effect on g can be substantial because of its exponential dependence and, in 
principle, it can easily overwhelm the geometrical piezoresistive effect. 

Of course, this tunneling mechanism of strain dependence should hold true also for other granular metals than 
the thick- film resistors as long as tunneling is an important element in transport properties. However the magnitude 
of the intrinsic piezoresistive effect depends on several factors such as elastic heterogeneity, diffusion of the metallic 
phase into the glass etc.. 

In this paper, we point out that a factor affecting in an important way the piezoresistive response is the degree of 
tortuosity the current has in flowing through a sample of granular metal. As an example, consider Fig. 1 where we 
report the results of a Monte Carlo calculation of a two dimensional bond-percolation resistor-network model. In the 
network, a fraction p of bonds distributed at random has a finite conductance g, while the remaining fraction 1 ~ p 
has zero conductance. When all bonds of the network are conducting {p — 1, Fig. la), the current flows along the 
direction of the applied voltage difference, that is the x direction. For p < 1 the missing bonds force the current to 
aquire nonzero components also along the direction perpendicular to the applied field (the y direction) as shown in 
Fig. lb where p = 0.8. For values of p close to the percolation thresold pc {pc = 0.5 for a simple square lattice), the 
tortuosity of the current flow is so developed that the contributions of components parallel and perpendicular to the 
applied voltage drop become comparable. This is clearly shown in Fig. Ic where the current flow is calculated for 
p = 0.55. 

How the tortuosity of the current flow affects the piezoresistive response, or better its anisotropy, is qualitatively 
described as follows. Immagine that each single bond in the directipa x of the applied field is slightly stretched so 
that its conductance becomes gx = 9 + g'x^, where e is the strain.O For simplicity, we keep the bonds along the 
direction perpendicular to the field unchanged: gy = g. The variation of the total conductance in this situation leads 
to the longitudinal piezoresistive response (applied strain parallel to the applied field). Such longitudinal response is 
nonzero for all the three cases depicted in Fig. 1 since there are always nonzero components of the current along the 
direction of the applied strain. Consider now the transversal piezoresistive response for which the strain is applied 
along the direction y perpendicular to the field, i.e., gy = g + g'yS and gx — g- Contrary to the longitudinal case, for 
p = 1 there are no components of the current along the y direction and the total conductance is therefore unaffected 
by the applied strain. In this case the transersal piezoresistive response is zero. However, for p < 1 current starts 
to flow also along the y direction and becomes affected by the strain induced variation of gy leading to a nonzero 
piezoristance effect. Since close to the percolation thresold the components of the current along the x and y directions 
are comparable (Fig. Ic), it is expected that in this case the longitudinal and the transversal piezoresistive responses 
have similar values. 

From the above qualitative discussion, we argue that away from the percolation thresold the piezoresistive response is 
highly anisotropic with the longitudinal component much larger than the transversal one, while close to the percolative 
critical point the two components becomes comparable, leading to an isotropic piezoresistive response. By using Monte 
Carlo calculations and results from the effective medium theory we show in the following that for resistor network 
models this conclusion holds true and that actually the longitudinal and the transversal piezoresistive responses are 
equal at the percolation thresold. Hence, a measurement of the piezoresistive anisotropy could be an alternative way 
to determine the critical concentration Xc of the conducting phase in granular metals. 

The paper is organized as follows. In the next section we introduce the model and, by using the effective medium 



theory, provide analytical formulas for the longitudinal and transversal piezoresistive responses. In Sec. HI we report 
our Monte Carlo calculations and compare them with the analytical formulas of the effective medium theory. The 
last section is devoted to a discussion and to the conclusions. 



II. EFFECTIVE MEDIUM THEORY 



In this section, we evaluate the piezoresistive response of a resistor network model by using a generalized effec- 
tive medium theory. Before entering the details of our model, let us first introduce an operative definition of the 
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piezoresistive effect and the different gauge factors. In Fig. 2 we show a schematic apparatus for the measurement of 
the piezoresistive response. The sample is placed on top of a cantilever bar clamped at one end, and a force F acts 
downward on the opposite end. We assume that in the absence of the force, transport is isotropic with conductance 
G. When F ^ 0, a tensile stress is built on the upper surface of the cantilever and a strain Exx = s is transfered to 
the sample. This strain is directed along the main axis of the cantilever which we choose to be directed along the 
X direction. If the thickness of the cantilever is sufficiently small compared to its width, the strain along y can in 
first approximation be neglected. If we assume that the strain field is completely transfered to the sample (which is 
a good approximation as long as the linear dimensions of the sample are much smaller than those of the cantilever) , 
then the strain field acting on the granular metal is: Sxx = £, £yy = 0, and Ezz = ~ i^), where v is the Poisson 

ratio of the material. By referring to Fig. 2, the longitudinal piezoresistive response is obtained by setting a potential 
difference between points a and b, and the total conductance up to the linear term in the strain is: 

Gx^G + &Gx^G{l~eGYi^), (4) 

where we have defined the longitudinal gauge factor: 

When the potential difference is applied between points c and d of Fig. 2, that is perpendicular to the strain in the 
X direction, the conductance can be expressed as Gy = G(l — eGFt), where GFt is the transversal gauge factor: 

Although GFl and GFt are the most used parameters, there is also a vertical gauge factor obtained applying a 
voltage difference between the upper and lower surfaces of the sample: GFz = —SGzjsG. 

Let us describe now our resistor network model. We consider a three dimensional simple cubic network where 
a bond has either a finite conductance g with probability p or zero conductance with probability 1 — p. We limit 
ourselves to the case in which the activation energy is sufficiently low and the temperature sufficiently high to 
approximate by a tunneling exponential. Furthermore, we allow a certain distribution in the tunneling distance in 
such a way that g — go exp[—K— AK{p' — 1/2)], where go is a prefactor, p' is a random variable uniformly distributed 
between and 1, and AK can assume values between zero (single tunneling distance) and AK = 2K. For this latter 
case therefore the tunneling exponent is continuously distributed between and -f^max = '^K. 

The probability distribution p{g) of the conductance value of a bond is: 
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Pig) = A^^h" -9] e^g-goe-«-""/^J +(l-p)%), (7) 

where Q{x) is the Heaviside step function. For a simple cubic network, the effective medium equation reads:0 

dgpig)^^ = 0, (8) 
9 + ^9 

which, by using Eq.(^, leads to the following expression for the effective macroscopic conductance g: 

As expected by the effective medium approximation, g vanishes for p — pc — 1/3, which could be interpreted as the 
effective field metal-insulator transition point. 

Now, Let us consider how the effective medium theory should be generalized to describe the effects of applied strain 
fields. To simplify the problem, we make the assumption that the strains are oriented along the directions of the 
bonds in the simple cubic network. Hence, in the presence of the strain field, the conductances become gx^ gy, and 
gz according to whether their corresponding conducting bonds are oriented towards to the x, y, and z directions, 
respectively. The probability distribution pi{g) of a bond in the i direction {i = x, y, z) is then equal to Eq.(^ where 
K and AK are replaced by Ki and AKi, respectively. Therefore, if en = 5di/di is the strain along the i direction, 
Ki = K{1 + £ii) and AKi — AK{1 + en). Within this model we have made the assumption that the strain field does 
not affect sensibly the prefactor go and from now on we set it equal to unity. The strain-induced|-a|nisotropy leads to 
a system of three coupled effective medium equations for the three macroscopic conductances gi :t2l 
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(10) 



where, by following the work of Bernasconi,Ej Si = 1/Ri — gi and Ri is the total resistance between two neighboring 
nodes in i direction for the effective lattice. For the x component, Rx reads: 



dx' dy' dz' ^ - 

- cos(a;')] + gy[l - cos(y')] + 9z[^ - cos(z')] ' 



(11) 



and Ry and Rz are obtained from ( pi]) by replacing in the numerator cos(a;') with cos(y') and cos(z'), respectively. 
For general values of the gi, the integrals Ri can be evaluated only numerically.tS However, we are interested in the 
piezoresistive effect which is just the linear response to an applied strain field so that, to our purposes, it is sufficient 
to evaluate the integrals Ri at the first order in the strain. As we show below this procedure permits an analytical 
solution of the anisotropic effective field equations. 

Let us replace in Eq.(ll) the effective conductances gi with = 5(1 ^ ^GFl), gy — 5(1 ^ ^GFt), and — 
g{l — eGFz), where, by following the definitions given above, we have identified the gauge factors for the effective 
conductances by: 



GFt 



£5 



GFt - 

£g 



GFz = - 



Sgz 
£g 



As shown in the Appendix, up to order e, R^ reduces to: 



Rx = 
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where ^ : 
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0.5264. Similarly, we obtain for Ry and R^: 
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(12) 

(13) 

(14) 
(15) 



By substituting Eqs.( |l3| - |l5| ) into Eq.(|lO|), and equating to zero the term linear in e, we obtain three linear equations 
for the three gauge factors. After some algebra, we obtain for the longitudinal and transversal components: 
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[3pAKX - 6KpY - AK][{1 - 2iy)^{9pY - AK) - l2pY{l - ly)] 
18j3r(l - ;/)[(4 - 9^)pY + ^AK] ' 



(16) 



GFq 



^[9pY - AK] [3pAKX - 6KpY - AK] (1 - 2iy) 
18py(l ~ t/)[(4 - 90pY + ^AK] ' 



where: 
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We show in the next section how GFl and GFt behave for different values of p, K, and AK in comparison with 
Monte Carlo numerical calculations. Before doing so, we stress here tVip main result of the effective medium theory, 
namely, that GFl and GFt are generally different unless p = Pc = l/3.Efl In fact, from Eqs.(p^p^ we obtain: 



lim GFl 
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(20) 



The effective medium theory therefore confirms the qualitative arguments discussed in the introduction and, as we 
shall see in the next section, it captures the essential physics governing piezoresistive anisotropy. It is also worth 
stressing that, at p = 1/3, GFl — GFt independently of the particular strain field used. 
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III. NUMERICAL CALCULATIONS 



In this section we compare the resuhs of the effective medium theory, Eqs.(|l6|-[19D, with numerical calculations 
of a three dimensional random resistor network model. In our Monte Carlo calculations, we have considered cubic 
networks with number of nodes up to 30 x 30 x 30 and bond conductances as described in the last section. The 
total conductance has been obtained by solving numerically the Kirchoff equations for all nodes in the network, and 
periodic boundary conditions have been imposed to the sides of the network where no voltage difference is applied. 
The different gauge factors have been obtained by calculating, for a fixed configuration of resistors, the difference in 
conductances when e = and e — 0.001. The results have been then averaged over a minimum of 30 to a maximum 
of 100 different runs. 

In Fig. 3 we show the results of Monte Carlo calculations (symbols) and effective medium theory (lines) for the 
limiting case for which there is no distribution in tunneling probabilities [AK = 0). In this limit, the effective medium 
formulas (p^p^ simplify considerably: 

hm GF, . ^3m-ri(1^2.) 4(1 3p)(l-.) 

lim GF, ^ K, (22) 

AA'^o ^ (l-i^)[9C(l-p)-4(l-3p)] ^ ' 

The results of Fig. 3 confirm the qualitative analysis made in the introduction (see also Fig. 1). Namely, for p = 1 the 
current flows exclusively along paths parallel to the direction x of the applied field leading to a non zero longitudinal 
piezoresistive effect (GFl = K) and a vanishing transversal response (GFt = 0). For p < 1, the missing bonds 
force the current to flow also along directions perpendicular to x. In this case, the piezoresistive response acquires 
a transversal component while the longitudinal one is lowered. Note that the analytical and numerical results agree 
quite well for bond probabilities larger than p ^ 0.6 while for lower values of p the effective medium results start to 
deviate from the Monte Carlo data. This is of course due the increased inaccuracy of the effective medium theory as. 
the percolation thresold is approached. However, at the critical bond probability {pc — 0.25 for simple cubic lattice,t3 
and pc — 1/3 for the corrsponding effective medium approximation) both the Monte Carlo data and the analytical 
results predict equal longitudinal and transversal responses. Finally, note that for v = 0.5, GFt remains equal to 
zero also for p < 1. This is due to the fact that, for this value of the Poisson ratio and our assumption on the applied 
strain field, one has Szz = — e, therefore the bonds in the z direction are contracted of the same amount as the bonds 
in the x direction are stretched. 

The results of Fig. 3 show clearly how piezoresistive isotropy is related to the vicinity to the percolation thresold. 
Another limiting situation is the one depicted in Fig. 4 where all bonds are present (p — 1) but with a finite ditribution 
of tunneling probability. Here, AK is varied between AK — and AK — 2K, K has been fixed equal to 4 and = 0, 
0.3, and 0.5. For AK = we recover the limit p = 1 of Fig. 3, while for finite values of AK the distribution of 
different tunneling probabilities induces tortuosity in the current flow which is reflected in a reduction of GFl and an 
enhancement of GFt. As before, the limit = 0.5 has the peculiarity of having GFt = independently of the current 
tortuosity. Note that, since for K = 4 and p = 1 the system is far away from the percolation thresold, effective medium 
theory agrees excellently well with the Monte Carlo calculations. This agreement becomes however less satisfactory 
as the distribution width increases as shown in Fig. 5 where results are reported for AK = 0, K, 2K as a function K 
for p = 1 . High values of tunneling distribution widths increase the fluctuactions of bond conductances reducing the 
validity of the effective medium approximation. 

Having analyzed two representative cases (p < 1 with AK — 0, and p = I with AK > 0) we show in Fig. 6 
the results for both p < 1 and AK > with K = A and i/ = 0.3. The difference between the longitudinal and 
transversal piezoresistive effects decreases as the bond conductance fluctuations increases whatever is the origin of 
such fluctuations (missing bonds or finite distribution of bond conductances). Note that, contrary to the other cases, 
for AK = 8 GFl increases as p decreases. From the effective medium formulas it is found that this feature persists 
also for higher K values if AK is somewhat larger than K. However it is difhcult to test this tendency with our 
Monte Carlo calculations. since we have employed the relaxation method for which large values of AK seriously slow 
the convergency of the recursive algorithm. 

The effect of bond fluctuations on the piezoresistive anisotropy is more clearly shown in Fig. 7 where we plot, for 
the Monte Carlo data of Fig. 6, the quantity x defined as: 

which measures the degree of anisotropy of the piezoresistive response. For a fixed value of AK, x decreases as p —> Pc 
from above, and it is expected to vanish at p = pc- x decreases also at fixed p as AK is made larger. We believe that 
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X, because of its high-p sensitivity for p ^ pc, could be a useful practical way to measure the proximity of a granular 
metal to its percolation critical point. 



IV. DISCUSSION AND CONCLUSIONS 



The results shown in the previous section clearly indicate how the anisotropy of the piezoresistive effect depends 
on the vicinity to the percolation thresold. We have interpreted the reduction of anisotropy as due to the increasing 
of bond conductance fluctuations as p — > Pc from above. Although our results are quite general, their application to 
real granular metals needs some additional comments. 

In our bond percolation model, we have assumed that the external strain fields are applied along the directions of the 
bonds in the network. However, in a more realistic situation the bond directions should be considered as random. In a 
practical calculation this could be achieved by considering an ensemble of ordered networks with different orientations 
with respect to the applied strain fields. The average over such an ensemble should be a satisfying description of 
a realistic case. What we expect, and bond-average effective medium results (not reported here) confirm, is that 
even for AK = and p = 1, the anisotropy parameter x, Eg. ([2^), is less than one. This is due to the random 
bond orientation which contribute, even if the network have strictly equal bond conductances, to a lowering of the 
piezoresistive anisotropy. However, also in this case, the reduction of x as p is lowered (or AK is enhanced) and the 
limit X = for p Pc remain valid. 

Another effect which has been ignored in the present analysis is the possibility of having elastic heterogeneity within 
the granular metal. A large difference between the elastic constants of the metallic and insulating phases can give 
rise to important fluctuations in the local strain fields. Since the microscopic tunneling processes are affected by the 
local rather than the macroscopical strain values, the elastic heterogeneity can influence the piezoresistive response 
in an important way. This is actually what is expected in thick-film resistors where very stiff metallic granules are 
embedded in a relatively elastically soft insulating glass. This large elastic heterogeneity could be at the origin of the 
large values of piezoresistive responses actually observed for this class of granular materials.cJ We expect however that 
the inclusion of this effect should infiuence the absolute values of GFl and GFt but leave the anisotropy parameter 
X relatively unaffected. 



APPENDIX A: 



In this appendix we report the evaluation of the integrals Rx, Ry, and Rz up to the linear term in the strain e. 
Let us consider first R^, Eq.(ll), and substitute the effective conductances gi, i = x,y,z, with = g{l — eGFl), 
gy — g{l — eCFx), and gz = g{l — eGFz), where the gauge factors are defined in Eq. (p!^) . Hence, up to the term 



linear in e we obtain: 



where 



Rx = ^ + ^GFi^e+UGFT + GFz)e, (Al) 
9 9 9 



1 r r r . . . 1 - cos(a;') 

n\.' , (A3) 



^0 ^0 



[3 — cos(a;') — cos(j/') — cos(z')] 



71-^ Jo Jo Jo [3 - cos(a;') - cos(y') - cos(z')]^ 

It is easy to show that 6 + 2c = a and a = 1/3. Moreover b can be rewritten as b — a — where: 

1 />7r pTZ />7r -1 

Form Ref . we have: 



^^^Jo Jo Jo 3-cos(x0-cos(,')-cos(z0- ^^'^ 



^ = A[i8 + 12V2 - lOVs - 7\/6]K[(2 - \/3)(\/3 - V2)f ~ 0.5264, (A6) 
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where K is the complete elhptic i nteg ral of the first kind. Equation (|l^) is obtained by substituting a = 1/3, 



b={l- ^)/3, and c = ^/6 into Eq.(|AlD. 
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FIG. 1. Calculated bond currents for a two-dimensional bond percolation model. The external voltage difference is applied 
between the left and right end sides of the network. The bond probabilities are p = 1 (a), p = 0.8 (b), and p = 0.55 (c). The 
directions and sizes of the arrows represent the directions and the intensities of the bond currents, respectively. 



FIG. 2. Schematical device for the measurement of the piezorosistivc response. The cantilever beam is clamped at one end 
and a force F acts downwards at the opposite end. The resulting strain field is transfered to the material of interest (dark 
region) . The longitudinal response is obtained when a potential difference is applied between points a and b and the resulting 
conductance is measured iov F ^ Q and F = 0. Instead, the transversal response is obtained when the potential difference is 
applied between c and d. 
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FIG. 3. Longitudinal and tranversal gauge factors as a function of bond probability p for different values of the Poisson ratio 

V. The tunneling probability has zero distribution {A.K — 0) and K = 4. Symbols refer to Monte carlo calculations (simple 
cubic lattice of 30 x 30 x 30 sites) while the solid and dashed lines are the results of the effective medium approximation. 
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FIG. 4. Longitudinal and tranversal gauge factors as a function of the tunneling distribution width A/sT for if = 4 and p = 1. 
Prom top to bottom: z/ = 0.0, 0.3, and 0.5. Symbols and lines have the same meaning as in Fig. 3 
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FIG. 5. Longitudinal and tranvcrsal gauge factors as a function of K for p = 1, v = 0.3 and different values of the tunneling 
distribution width Aif . Symbols and lines have the same meaning as in Fig. 3 



12 




0.0 0.2 0.4 pO.6 0.8 1.0 



FIG. 6. Longitudinal and tranversal gauge factors as a function of p for different values of the tunneling distribution width 
AK. K = A and v = 0.3. Symbols and lines have the same meaning as in Fig. 3 
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FIG. 7. Factor of piezoresistive anisotropy x, Eq.(|23|), for the Monte Carlo results of Fig. 6. 
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